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Abstract 

Smooth double crossing of the phantom barrier uja = — 1 has been found possible in cosmological model 
with Gauss-Bonnet-scalar interaction, in the presence of background cold dark matter. Such crossing has 
been observed to be a sufficiently late time phenomena and independent of the sign of Gauss-Bonnet-scalar 
interaction. The luminosity distance versus redshift curve shows a perfect fit with the ACDM model up to 
z = 3.5. 
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c3 1 Introduction 

] The puzzle associated with recent cosmic acceleration, triggered by 70% of dark energy or more Ij is far from 
^ ■ being resolved uniquely. In the mean time, cosmologists are being confronted with yet another more intriguing 
' challenge to explain the crossing of the so called phantom divide line {wa = —1), at sufficiently late time of 
l/^ cosmological evolution. Some recent analysis [I], [2] of the presently available observational data are in favour of 
the value Wde < — 1, at present., Wde being the dark energy equation of state. There are also a lot of evidence 
. . all around [3], of a dynamical dark energy equation of state, which has crossed the so called phantom divide line 
' Wa = — 1 recently, at the value of red-shift parameter z « 0.2. Apparently though the problem turns out to be 
more serious and complicated, but then, the puzzle of crossing the phantom divide line has also rendered some sort 
of selection rule. ACDM-model, which is known to suffer from the disease of fine tuning (see W for a comprehen- 
sive review) can now be ruled out due to the requirement of a dynamic state parameter. Further, if the analysis of 
Vikman [S] is correct, then it is not possible to cross the phantom divide line in a single minimally coupled scalar 
field theory, without violating the stability both at the classical [6^ and also at the quantum mechanical levels [7], 
?H (though it has recently been inferred [8] that quantum Effects which induce the w < —1 phase, are stable in the 
(p* model). Thus single minimally coupled scalar field models like quintessence (w > —1) and phantom {w < —1) 
are to be kept aside. Consequently, we are now left with some what more complicated models. One of these is a 
hybrid model, composed of two scalar fields, viz, quintessence and phantom - usually dubbed as quintom model 
[S]. Other models like non- minimal scalar tensor theory of gravity |llO , hessence PTj and models including higher 
order curvature invariant terms [12' also exist in the literature. 

Gauss-Bonnet term is yet another candidate which may be pursued for the purpose. The possibility of crossing the 
phantom divide line through Gauss-Bonnet interaction has been explored in some recent works [13] , [2] . But then, 
these models are even complicated in the sense that either brane-world scenario T3] or scalar field and matter 
coupling }14i are invoked. In this article the possibility of smooth crossing of the phantom divide line wa = — 1 has 
been expatiated simply by introducing Gauss-Bonnet-Scalar coupling term in the 4-dimensional Einstein-Hilbert 
action. 

Gauss-Bonnet term arises naturally as the leading order of the a' expansion of hcterotic supcrstring theory, 
where, a' is the inverse string tension [15]. Gauss-Bonnet term is topologically invariant and thus does not con- 
tribute to the field equations in four dimensions. However, the low energy limit of the string theory gives rise to 
the dilatonic scalar field which is found to be coupled with various curvature invariant terms |16| . The leading 
quadratic correction gives rise to Gauss-Bonnet term with a dilatonic coupling ;17J. Therefore it is reasonable to 
consider Gauss-Bonnet interaction in four dimension with dilatonic-scalar coupling. Several works with Gauss- 
Bonnct-dilatonic coupling are already present in the literature TS\. In particular, important issues like - late time 
dominance of dark energy after a scaling matter era and thus alleviating the coincidence problem, crossing the 



phantom divide line and compatibility with the observed spectrum of cosmic background radiation have also been 
addressed recently [19] . 

In a recent work with Gauss-Bonnet interaction [20] . a solution in the form a = age^^ (a being the scale 
factor, and A > 0) has is been found to satisfy the field equations with different forms (sum of exponentials, 
sum of inverse exponentials, sum of powers and even quadratic) of potentials. Solution in a more general form 
(a = aoe^* ),A > 0,0 < f < 1, for Einstein's gravity with a minimally coupled scalar field was found in the 
nineties [21] and was dubbed as intermediate inflation. We [20|, on the other hand, observed that such solution 
depicts a transition from decelerated to accelerated expansion at sufficiently later epoch of cosmic evolution, which 
asymptotically goes over to de-Sitter expansion. Thus, it appeared that such solution may construct viable cos- 
mological models of present interests. Under this consequence, a comprehensive analysis has been carried out [52] 
with such solution in the context of a generalized k-essence model. It has been observed that it admits scaling 
solution with a natural exit from it at a later epoch of cosmic evolution, leading to late time acceleration with 
asymptotic de-Sitter expansion. The corresponding scalar field has also been found to behave as a tracker field 
[23j , thus avoiding cosmic coincidence problem. 

In the present work, we show that Gauss-Bonnet-Dilatonic scalar coupling with Einstein's gravity in four dimen- 
sions, admits solution in a general form (a = aoe"**^), A > 0,0 < f < 1, which is viable of crossing the phantom 
divide line twice, once from above and the other from below in the recent epoch. Since the crossing is transient, so 
we may conclude that it does not show any pathological behaviour like Big- Rip [6], at least in the classical level. 

2 The Model with Gauss-Bonnet Interaction 

We start with the following action containing Gauss-Bonnet interaction 



/ 



d'^V^i^ + ^G{R) - i^,^ 0"^ - F(0) + L^], (1) 



where, 

GiR) - - AR^.R'"' + R^,,p„R'"'P" 

is the Gauss-Bonnet term which appears in the action with a coupling parameter A((/)) and Lm is the matter 
Lagrangian. For the spatially flat Robertson- Walker space-time 

ds^ = ~dt^ + a^{t)[dr'^ + r'^de'^ + sin^ed<j)^], 

the field equations in terms of the Hubble parameter H = ^, are 

2H + 3H^ = -[^02 _ y(^) ^ 2K'^{HH + H^) + (A> + A"^'')H^ + p,n] = -{pde + Pm), (2) 

= [^02 + y(^) _ 3A>i/3 ^ ^ (^^^ + p^)^ (3) 

in the units k'^{= 8ttG) = H = c — 1 . In our analysis the Gauss-Bonnet scalar interaction plays the role of 
dark energy, for which suffix (de) has been introduced. Thus, pde and pde are the effective pressure and the 
energy density generated by the Gauss-Bonnet-scalar interaction, while Pm and pm are the pressure and the 
energy density corresponding to background matter distribution respectively. The background matter satisfies the 
equation of state, 

p™=P^a-3(i+-"), (4) 

where, pi is a constant and Wm is the state parameter of the background matter. In addition we have got the (f> 
variation equation 

{4> + 3iJ0 + V) = SA'H'^iH + H^), 



which is not an independent equation and will not be required in our analysis. In the above, over-dot and dash 
(/) stand for differentiations with respect to the proper time t and (j) respectively. Now, in view of equations (2) 
through (4), we are required to solve for a, 0, V{4>), A(0),p„i and p,„ , which requires three additional assumptions. 
Firstly, we consider that the Universe is filled with cold dark background matter with equation of state, pm = 
while the second assumption is the one made previously in [20], viz., 



A'0 = A, (5) 

where, A is a constant. This, as indicated in "20\ is physically reasonable, since it implies that the Gauss-Bonnet 
coupling parameter A{(j){t)) = Xt, grows in time to contribute at the later epoch of cosmological evolution. In 
view of the above assumption the field equations (2) through (4) are expressed as, 

2H + 3ij2 = - [1^2 „ v{(f,) + 2XHH + 2XH^] = -pde, (6) 

and, 

Pm = pia~^. (8) 
Now, for our third assumption, we start from the ansatz. 



with < / < l,n = A^^ > 0, which leads to the form of the solution of the scale factor mentioned in the 
introduction. Thus the complete set of solutions are given by, 



A . 3 f, 2/(1-/) 3/2 3/2 
a = aoexp( — ); p,„ = -3 exp(--i^ j; p^e = — 7(^^\ 57WT~n! Pde = y.^,. 



Wde = a] 



3 I 2n/(l - f)t-f - 3/2 



3alP - p,nH^i^-f) exp(-ft/) 



_ 6/(l-/) 6/2 P. ^ _ 2/(1 - /) p. 

Pde + .iPde- ^^(2_/) „2i2(i-/) alexpi^f)' P'^' + ' nt^~f a3exp(fi/)' ^ ' 

•2 A/3 , 2A/2(1-/) , 2/(1-/) p. 
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3A/3 A/2(l-/) , 3/2 /(I-/) 



2n3t3(i-/) n2t3-2/ „2i2(i-/) ^^2-/ 2a3exp(fi/)' 



Above set of solutions (10) indicates that such a model of the Universe admits an early deceleration, but during 
evolution it starts accelerating since strong energy condition is violated, pde + Spde < 0. Further, the dark energy 
equation of state also admits the possibility of crossing the wa = — 1 line, since, transient violation of the weak 
energy condition, pde + Pde < is seemingly possible. Finally, the equation of state asymptotically touches the 
wa = — 1 line from above and behaves as cosmological constant. To show such behavior graphically, let us express 
the state parameter Wde in terms of the red-shift parameter z which is defined as. 



l + z = %^=exp[i(t/-t/)], 



a{t) 




z 



Figure 1: State parameters Wdeiz) has been plotted against the red-shift parameter z, (with, / = 0.5, ft. = 
0.66, to = 13 Gyr,i}jno = 0.26). Smooth double crossing of the Cosmological constant barrier is observed at 
sufficiently later epoch, z w 1.92 from above and z « 0.39 from below. 



where, a{to) is the present value of the scale factor, while a{t) is that value at some arbitrary time t, when the 
light was emitted from a cosmological source. Thus, 

= t{ - 7iln(l + z). (11) 
In view of equation (11) Wde can be expressed as, 



{2nf{l- f)~3f[tf-n\nil + z)] 
^ 3/2 [tl - n In (1 + z)] - p,n^ [tl - 7i In (1 + z)] ^ exp {- ^[tl - nln {1 + z)]) ^ 



^de = I —7-7 — 7 — — -2^7 —m — — I ' (12) 



where, Uq has been set equal to one without any loss of generality. Let us now choose / — 0.5. The motivation of 
choosing such a value of / is twofold. Primarily, it is impossible to find an explicit form of the potential V = V{cf)) , 
otherwise. Further, since n has the dimension of t-^ , so the parameter gets a comfortable dimension of time. 
If we now take up some more numbers, like the present value of the Hubble parameter and the age of the 

Universe to as, 

H-'^^Gyr, t„ = 13 Gyr, 

then, for h = 0.66, n can be found from the ansatz (9) a,s n — 2.0552. Further taking the present value of the 
matter density parameter flmo = 0.26, we find in view of solution (10), 

o P"'" , go~^exp(-ft;^) 

"mo = = Pi{ 7; ) = 0.26, 

Pco 3 

where, pmo and Pco are the present values of the matter density and the critical density respectively. Thus, we 
find, 

n^Pi = 2.897. 

Noting that in this model the red-shift parameter does not go beyond the value z = 4.78, we plot the dark energy 
equation of state parameter (wde) versus the red-shift parameter (z) in figure (1). It is apparent that the phantom 
divide line lua has been crossed twice, once from above at z « 1.92 and then from below at z « 0.39. Such double 
crossing of the phantom divide line is devoid of any sort of pathological behaviour. 

To check how far our present model fits with the standard KCDM model, we also make the luminosity-redshift 
and distance modulus-redshift plots. For KG DM model the relation between the luminosity distance and redshift 
is the following, 

dz 

H,dL = {l + z)j^ 0y,^0 26(l + z)3 ' 




Figure 2: The fit is almost perfect up to z = 3.5. A little discrepancy is observed there after as KCDM model 
(blue) slightly takes over the present model (red). 




Figure 3: The fit is absolutely perfect and the two models (blue for KCDM and red for the present) are practically 
indistinguishable. 



while the relation corresponding to the present model is, 
1 + z 

HqcIL = — ^— / [\/tQ ~ nln {\ + z)]dz , 
Vto Jo 

with to = 13 Gyr., and n = 2.055. The plot (figure 2) shows a perfect fit between the two models up to z = 3.5. 
There is a little discrepancy there after. Since, luminosity distance has already been expressed as a function of 
redshift, so the relation between distance modulus and redshift may be found in view of the following equation, 

m-M = 51ogio(^) + 25, 

where, m and M are the apparent and absolute bolometric magnitudes respectively. However, since we use HqcII 
instead, so our relation is slightly modified as, 

m- M ^ 51ogio(DL) + 31, 

where, DL — H^dL . The plot (figure 3) demonstrates that the two models are practically indistinguishable. 

Now we can proceed to make some even more comfortable choice of the parameters of the theory, like / = 0.5 , 
and to = 13 Gyr. as before, but with n = 2, for which = 14.42 Gyr., in view of ansatz (9), which corresponds 
to /i w 0.68. The scale factor now has a convenient form as a = exp . With these values one can find. 



n^p, = 3.09, 




Figure 4: State parameters Wdeiz) has been plotted against the red-shift parameter z, (with, / = 0.5, ft. ~ 
0.68, to = 13 Gyr.,flmo — 0.24). Smooth double crossing of the Cosmological constant barrier is observed at 
sufficiently later epoch, z w 2.085 from above and z w 0.46 from below. 
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Figure 5: The fit is perfect up-to z = 3.6. The ACDM model (blue) slightly overtakes the present model (red) 
there after. 



for ri„jo — 0.24. The plot (fig. 4) is almost the same as before with transient double crossing. Other plots viz., 
luminosity distance versus redshift (in figure 5) and distance modulus versus redshift (in figure 6) show even better 
fit than the earlier one, with the standard ACDM model. 

Thus, we observe that with the age to — 13 Gyr., and 0.66 < h < 0.68, such transient crossing of the phantom 
divide line is permissible for the present value of dark energy density parameter ^de\present > 0.74. In order to 
consider some higher value of the age of the Universe, to = 13.73 Gyr. (say) as suggested by Spergel et alp4]. 
either one has to go to almost the lowest limiting value of ft, « 0.61 [25] or one has to accept much higher value of 
the present dark energy density parameter ^de\present > 0.78, otherwise, the state parameter versus redshift plot 
shows certain discontinuities. We certainly remember that in order to simplify the field equations considerably, we 
have made one important assumption, viz., in equation (5). Relaxing this assumption one might get rid of such 
discontinuities, as well. This we pose to study in a future communication. 

So far, we remain silent about the form of the potential. It is simply because, despite the most convenient choice 
of the parameter, / = 0.5, it is still impossible to find an analytical solution for (f>, in view of the solution (10). 
As a result, the form of the potential as a function of (j) remains obscure. However, we can plot the potential 
as a function of time, by choosing our second case, n — 2, for further simplification. It is important to note 
that though the results sofar obtained, are independent of the value and signature of A , the form of the potential 
depends largely on it. In the following we make three such plots (taking the help of "Manipulation" programme 
of Mathematica 6) to show how the form of the potential changes with different values (starting from negative to 
large positive) of A . 
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Figure 6: The fit is absolutely perfect and there is practically no way to distinguish the KCDM model (blue) 
with the present one (red). 




Figure 7: The form of the potential as a function of time for A < —9 . Note that it remains negative throughout 
the evolution. 




Figure 8: The form of the potential for A w 65 
future. 



Note that it is zero at the present epoch but tends to grow in the 
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Figure 9: The form of the potential for A = 200. Note that the form is appreciably different. 



3 Concluding remarks 

Altogether we have obtained a late time transient crossing of the phantom divide line first from above and more 
recently from below, starting from the inclusion of a Gauss-Bonnet-dilatonic scalar coupling term in the standard 
Einstein-Hilbert action in four dimension. Since the crossing is transient, so such double crossing is free from any 
sort of pathological behaviour both at the classical [5] and at the quantum mechanical levels. The striking 
feature of the model lyes in it's indistinguishability with the standard AC DM model, in terms of the luminosity- 
redshift and more precisely for distance modulus-redshift curves. To identify between the two models we therefore 
require to observe dark energy equation of state Wde independently. If Wde is truly found to be dynamical and 
has really encountered a recent crossing of the phantom divide line, then only we can definitely distinguish the 
standard ACDM model with the present one. It is highly interesting to learn that smooth transient crossing of 
the phantom divide line is allowed for both negative and positive (A ^ 0) type of Gauss-Bonnet-scalar interaction. 
Figures (7), (8) and (9) also reveal that it is true even for different forms of the potential. Such transient crossing 
independent of the value of A also signals that it might be possible to carry out the same treatment even for a 
single scalar field model. 
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